Factorization in Break-up and Recombination Processes 
for Atoms with a Large Scattering Length 
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Break-up and recombination processes for loosely-bound molecules composed of atoms with a 
large scattering length a necessarily involve interactions that are nonperturbative in the exact 2- 
body interaction. If these processes involve atoms with relative momenta much larger than 7i/|a|, the 
leading contributions to their rates can be separated into short-distance factors that are insensitive 
to o and long-distance factors that are insensitive to the range of the interaction. These factorization 
contributions can be obtained from the leading term in a perturbation expansion in the exact atom- 
atom scattering amplitude. The short-distance factors are atom-atom cross sections at a lower 
collision energy. In the special case of inclusive break-up cross sections for atom-molecule scattering, 
the long-distance factors simply count the number of atoms in the molecule. 

PACS numbers: 12.38.-t, 12.38.Bx, 13.20.Gd, 14.40.Gx 
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Nonrelativistic particles with short-range interactions 
that have been tuned, either by the experimenter or 
fortuitously by nature, so that their S-wave scattering 
length a is much larger than the range, have universal 
low-energy properties that depend on a but are other- 
wise insensitive to their interactions at short distances. 
(See Ref. pj and references therein.) If the particles form 
loosely-bound 2-body or higher TV-body clusters whose 
sizes are comparable to \a\, the clusters also have uni- 
versal properties. A classic example in atomic physics is 
4 He atoms whose scattering length a » 100 A is much 
larger than their effective range r s 7 A. The 4 He dimer 
and the excited state of the 4 He trimer both have sizes 
comparable to a. A classic example in nuclear physics is 
nucleons. The deuteron is a bound state of the neutron 
and proton associated with a large scattering length in 
the spin-triplet channel. 

The large scattering length implies that interactions 
between atoms whose relative momenta are comparable 
to h/\a\ are nonperturbative in the 2-body scattering am- 
plitude. A simple consequence of these strong interac- 
tions with a > is the existence of a loosely-bound dimer 
whose binding energy is given by 



h 2 /(ma 2 ). 



(1) 



In some cases (for example, identical bosons), the strong 
interactions produce the Efimov effect: as a — ► ±oo, there 
are increasingly many loosely-bound trimers with an ac- 
cumulation point at the 3-atom threshold . In the limit 
a = ±oo, the trimers have an asymptotically exponential 
spectrum: 
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where sq is a numerical constant whose value for iden- 
tical bosons is 1.00624, rt* is an arbitrary integer, and 
is a 3-body parameter pj. The strong interactions 
can also lead to intricate dependence on the scattering 



length a. An example in the case of identical bosons is 
the event rate constant in the low-energy limit for 3-body 
recombination into the loosely-bound dimer: 



A : 



7687r 2 (47r - 3V3>a 4 /7 



smb. (ttsq) + cosh (ttsq) cot 2 [so ln(a/t*) + 1.16] 



(3) 

The log-periodic dependence on a in Eq. (PJ) was discov- 
ered in Refs. |3j. The completely analytic expression was 
derived more recently by Petrov Q. The coefficient of 
ha 4 /m on the right side of Eq. © can range from to 
402.7 depending on the value of the 3-body parameter k*. 
The dimer break-up cross section in atom-dimer scatter- 
ing at a collision energy E just above the threshold Ery 
has the same log-periodic factor: 
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To reproduce such intricate dependence on a requires 
accurate numerical methods that are capable of solv- 
ing the problem "exactly", that is, to any desired pre- 
cision. The Schrodinger equation for 2 or 3 atoms in- 
teracting through a short-range potential can be solved 
exactly. Systems of 4 atoms are at the frontiers of few- 
body physics: the binding energies of tetramers can be 
calculated accurately, but there are some scattering ob- 
servables for which effective calculational methods have 
yet to be developed. In any process that involves loosely- 
bound molecules, strong interactions involving the mo- 
mentum scale h/\a\ are unavoidable. This suggests that 
accurate calculations of the rates for such processes may 
be possible using present methods only if the total num- 
ber of atoms involved is at most 3 or 4. 

In this Letter, we point out that this obstacle can be 
avoided for processes involving atoms whose momenta 
relative to the loosely-bound molecules are set by a scale 
Q 7i/|a|. Examples of such processes are the break- 
up of the molecule in a collision with an energetic atom 
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and the formation of the molecule in a collision involv- 
ing energetic atoms. The separation of scales Q 3> h/\a\ 
can be exploited by using factorization to separate the 
rate into short- distance factors that are insensitive to a 
and long-distance factors that do not involve the scale 
Q. The strong interactions associated with large a enter 
only in the long-distance factors. The rate can be cal- 
culated accurately using present methods provided the 
long-distance factors involve at most 3 or 4 atoms. In 
some cases, even this limitation is unnecessary because 
the long-distance factors can be determined analytically. 
Thus factorization significantly expands the list of pro- 
cesses whose rates can be calculated accurately using 
present methods. 

Factorization has proved to be a powerful tool in quan- 
tum chromodynamics (QCD). The Nobel Prize in physics 
for 2004 was awarded to Gross, Politzer, and Wilczek 
for the discovery of the asymptotic freedom of QCD ||. 
Asymptotic freedom refers to the decrease of the strength 
of the interaction between quarks as their separation de- 
creases, or equivalently, as their relative momentum in- 
creases. For selected observables, asymptotic freedom 
can be exploited by using factorization to separate the 
observables into short-distance factors that involve only 
weak QCD interactions between quarks and gluons and 
long-distance factors that involve strong QCD interac- 
tions [6|. The short-distance factors involve a large mo- 
mentum transfer Q and can be calculated using pertur- 
bation theory in the running coupling constant a s (Q) 
of QCD. The long-distance factors involve only momen- 
tum transfers small compared to Q. In most cases, sys- 
tematic methods for calculating the long-distance factors 
have not been developed. The bulk of the quantitative 
evidence that QCD describes the strong force is the ex- 
perimental verification that the short-distance factors are 
correctly predicted by perturbative QCD. 

Atoms with a large scattering length have a property 
analogous to asymptotic freedom: the strength of their 
interaction, as measured by the magnitude of the elastic 
scattering amplitude, decreases as their relative momen- 
tum hk increases. This behavior of the atom-atom in- 
teraction motivates our use of the factorization strategy 
that has been so successful in QCD. The S-wave elastic 
scattering amplitude has the form 



f k = [(-l/a+±r s fc 2 + ...)-ik] 1 



(5) 



where r s is the effective range. For collision energies E at 
which higher partial waves can be neglected, the elastic 
cross section for identical bosons is 
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The cross section has its maximum value 8na at E = 0. 
In the scaling region h 2 /mo? <£< Ti 2 /mr 2 , the cross 
section scales like 1/E: 
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FIG. 1: Feynman diagram for the break-up of a dimer by the 
hard scattering of one of its constituent atoms and a colliding 
atom. The blob represents the wavefunction of the dimer and 
the dot represents the exact atom-atom scattering amplitude. 



Consider the break-up of a loosely-bound dimer 
through a collision with an atom. We take all three 
atoms to be identical bosons. The two atoms inside the 
dimer have typical separation a. The atom-dimer colli- 
sion energy is E — jQ 2 /m, where Q is the momentum 
of the atom or the dimer in the center-of-mass frame. If 
Q S> h/a, the leading contribution to the break-up of 
the dimer comes from one of its constituent atoms be- 
ing knocked out by the energetic incoming atom, while 
the other acts as a spectator. This process can be rep- 
resented by the Feynman diagram in Fig.^ which is the 
leading term in a perturbation expansion in the exact 
atom-atom scattering amplitude. The expression for the 
T-matrix element is 
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where En = h /{ma 2 ){\ + r s /a + . . .) is the dimer bind- 
ing energy including range corrections, qs is the relative 
momentum of the two atoms in the dimer (one of which 
is the outgoing atom labelled 3), and hkyi — |Q is the 
relative momentum of the incoming atom and the scat- 
tered atom from the dimer (which is also the relative 
momentum of the outgoing atoms labelled 1 and 2). The 
sum is over cyclic permutations of the three outgoing 
atoms. The effects of subsequent interactions between 
the scattered atoms and the spectator atom in the dimer 
enter through higher-order diagrams in the perturbation 
expansion in the exact atom-atom scattering amplitude. 
These diagrams are all suppressed by a factor of hj (aQ) 
and can be neglected if the collision energy is sufficiently 
high. If we take the limit Q ^> h/a, the T-matrix el- 
ement in Eq. JSJ can be expressed as the product of a 
short-distance factor that involves only the large momen- 
tum scales Q and h/r s and a long-distance factor that 
involves only the small momentum scales q and h/a: 
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The long-distance factor is proportional to the 
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momentum-space wavefunction of the dimer in the zero- 
range limit. The cross section is obtained by squaring 
the T-matrix element and integrating over the momenta 
of the three identical outgoing atoms. The cross sec- 
tion is dominated by the squares of the terms in Eq. JJjJ , 
which we call diagonal terms. The terms in the cross sec- 
tion that correspond to the interference between terms in 
Eq. @, which we call cross terms, are suppressed by a 
power of Ti/(aQ). We can take the limit Q 3> h/a in the 
momentum integrals. The result for the break-up cross 
section at collision energy E ^§> h 2 /ma 2 can be expressed 
in the factored form 
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The coefficient 2, which is the long-distance factor, is ob- 
tained by using the normalization integral for the dimer 
wavefunction. The factorization form in Eq. (|lCj|> applies 
at arbitrarily large E provided a aa is the total atom- 
atom cross section. In the scaling region Ed <C E <C 
h /mr 2 , this breakup cross section scales as 1/E: 
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The result can be generalized to the inclusive break-up 
cross sections for Efimov trimers and other loosely-bound 
molecules whose constituents all have separations of or- 
der \a\. The leading contribution comes from one of the 
atoms in the molecule being knocked out by a collision 
with the energetic incoming atom. If the molecule con- 
tains N loosely-bound atoms, the atom-molecule collision 
energy is E = Q 2 /m, where Q is the momentum in 
the center-of-mass frame. The short-distance factor in 
the T-matrix element is the same as in Eq. JJJJ, except 
that the relative momentum is hk = ^ff-Q- For any spe- 
cific final state, the long-distance factor may be a very 
complicated function of a and the 3-body parameter k* 
that appears in Eq. @- It involves the wavefunction of 
the molecule and the wavefunction of the (N — l)-atom 
system that remains after the hard scattering. Some of 
those N — 1 atoms may be bound into dimers or other 
clusters. However, if we sum over all these (N — l)-atom 
states, we can use completeness relations to show that the 
long-distance factor is simply the number N of atoms in 
the molecule. Our result for the inclusive break-up cross 
section as a function of the collision energy E is 



(break-up) 
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This expression is valid when the collision energy satis- 
fies E ^> h 2 /ma 2 and E ^> Em, where Em is the bind- 
ing energy of the molecule with respect to the iV-atom 
threshold. 

The cross section for the break-up of 4 He dimers in 
collisions with Xe atoms has been studied in Ref. |7j for 
collision energies E ranging from 46 K to 348 K. Since the 
Xe atom is so much heavier than a He atom, each of the 
atoms in the He dimer carries approximately ^ the col- 
lision energy E. The analog of our factorization formula 



in Eq. I|10f) predicts that if E is much larger than the He 
dimer binding energy, which is about 1.6 mK, the He2~Xe 
break-up cross section should be approximately 2 times 
the He-Xe cross section at collision energy \E. This 
prediction agrees reasonably well with the VCC-RIOS 
approximation Q studied in Ref. Q ■ By comparing the 
cross section cr^ C ' RIOS in Table 2 of Ref. with the 
appropriate interpolated values of cry cXc , we can see that 
the difference is less than 9% in the energy range from 93 
K to 348 K. In the Independent- Atoms model considered 
in Ref. Q, the total cross section is 2 times the He-Xe 
cross section at collision energy E instead of \E. This 
model underestimates the VCC-RIOS cross section by an 
amount that ranges from 11% at 93 K to 30% at 348 K. 

Factorization can also be applied to the rate for 3- 
body recombination into the loosely-bound dimer when 
the three incoming atoms all have relative momenta much 
greater than h/a. The leading contribution comes from 
a hard scattering of two of the atoms that scatters one 
of them into a state with small momentum relative to 
the third atom, followed by the coalescence of those two 
atoms into a dimer. The coalescence probability is sub- 
stantial only if the relative momentum of the scattered 
atom and the third atom is of order h/a. At leading 
order in the perturbation expansion in the exact atom- 
atom scattering amplitude, the T-matrix element for this 
process can be expressed as the sum of three Feynman di- 
agrams obtained by time-reversing the diagram in Fig.Q] 
and summing over cyclic permutations of the three in- 
coming atoms. We work in the center-of-mass frame, 
taking the momenta of the three incoming atoms to be 
Pi, P2, and f>3 and the momenta of the outgoing dimer 
and atom to be +Q and — Q. The collision energy is 
E = (Pi +P2 (2m) . For large collision energy E, the 

leading contribution to the T-matrix element is given by 
Eq. Each of the three terms in the sum is the prod- 
uct of a short-distance factor and a long-distance factor. 
The hard scattering of atoms 1 and 2 followed by the coa- 
lescence of one of them with atom 3 gives the term shown 
explicitly in Eq. where hki2 = \{pi — pi) is the large 
relative momentum of atoms 1 and 2 and <h = pz — \Q 
is the small relative momentum between atom 3 and the 
atom it coalesces with. The other two terms are ob- 
tained by cyclic permutations of 1, 2, and 3. The 3-body 
recombination rate is obtained by squaring the T-matrix 
element and integrating over the final momenta of the 
outgoing atom and dimer. In the limit E 3> h 2 /(ma 2 ), 
the rate R for forming dimers per volume and per time 
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6(E-3 P 2 /m). (13) 



The three terms in the sum are diagonal terms that come 
from the squares of the terms in Eq. The cross 

terms corresponding to the interference between terms in 
Eq. JjJ) have been neglected because they are suppressed 
by a factor of (Ed/E) 1 / 2 . 
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FIG. 2: Three-body recombination event rate K3 (in units 
of cm 6 /s) as a function of the collision energy E (in units of 
mK). The solid lines are the total result in Fig. 3 of Ref. [jj. 
The dashed lines are the hyperspherical averages of the R in 
Eq. I|170 . the diagonal terms, and the cross terms. The dotted 
line is the scaling approximation in Eq. (1161 . The horizontal 
arrow indicates the maximum possible value of Kz at E — 
allowed by Eq. 



The event rate R in Eq. (|13|) can be reduced to 
a function K3 of the collision energy E by averaging 
over the momentum hyperangle 0*3 denned by pk — 
(^mE) 1 / 2 cosafc and over the angle /3 3 between the vec- 
tors P3 and pi — pi. The hyperangular average of R can 
be expressed as 

(R) = - da 3 sin 2 (2a 3 ) / df3 3 sin (3 3 R. (14) 
* Jo Jo 

The hyperangular average of the rate in Eq. (|13|) is 

/D . 384\/37r?i 3 3 
K3 = (R) ~ j— <taa{\E). (15) 



This result is valid when E ^> En. In the scaling region 
E D < E < ?i 2 / mr 2 , the rate in Eq. (|15fl scales as 1/E 2 : 



12288V37T 2 ?i 5 
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Suno ct al. have calculated the 3-body recombina- 
tion rate of 4 He atoms into the dimer at collision en- 
ergies up to 10 mK using accurate solutions of the 3- 
body Schrodinger equation for the HFD-B3-FCI1 poten- 
tial model 0. In this model, the scattering length is 
a = 91.0 A, the effective range is r s ss 7 A, and the 
dimer binding energy is Ed = 1.6 mK. The scattering 
length deduced from the universal formula for the dimer 
binding energy in Eq. is an = 87.0 A. In Fig. 01 we 
compare the result of Ref. j£| for K 3 with the scaling 
approximation in Eq. (|16|) . At E — 10 mK, the highest 
collision energy considered in Ref. 0, our scaling ap- 
proximation in Eq. 1.1(1) is larger than the total result in 
Ref. [9j by a factor of 28.8. Our factorization approxima- 
tion in Eq. I|15fl is larger than the total result in Ref. @ 
by a factor of 23.7. 

We now consider why the scaling approximation in 
Eq. (|16fl overestimates the 3-body recombination rate at 
10 mK by more than an order of magnitude. We will 
show that the scaling approximation can be accurate only 
for collision energies that are at least an order of mag- 
nitude greater than 20 mK. We calculate the sum of the 
contributions to R at leading order in the perturbation 
expansion in the exact atom-atom scattering amplitude 
without making the factorization approximations. The 
T- matrix element is given by Eq. (JSJ. We set the effec- 
tive range r s to zero to make the scaling behavior of the 
individual contributions at large energy E more evident. 
The resulting expression for the 3-body recombination 
rate is 



m E 
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(cos 2 a 12 (t) - \ + \E D /Ef + \{E D /E){1 + E D /E) + 3t(l -<)(! + E D /E) sin 2 a 3 



(17) 



The angle ai 2 (i) is defined by 

cos 2 a\2 (t) = t cos 2 ai + (1 — t) cos 2 a 2 . 



(18) 



The two terms inside the sum in Eq. I|17|) are a diag- 
onal term, which corresponds to the square of a term 
in Eq. |JSJ, and a cross term, which corresponds to the 
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The hyperangular average of the rate R in Eq. (|17(l can 
then be calculated using Eq. (fH|) . In Fig. we show 
K3 = (R) as a function of the collision energy E, as well 
as the contributions to K% from the diagonal terms and 
from the cross terms. The cross terms are larger than 
the diagonal terms at low energies and smaller at higher 
energies, with the crossover occurring near E = 20 mK. 
At very high energy, the cross terms scale like E~ b l 2 . 
They eventually become negligible compared to the di- 
agonal terms, which scale like E~ 2 . The factorization 
approximation in Eq. 1)15(1 is a high energy approxima- 
tion to the contribution from the diagonal terms. From 
the crossover point and the scaling behavior, we can in- 
fer that the factorization approximation can be a good 
approximation only if the collision energy exceeds the 20 
mK by more than an order of magnitude. 



FIG. 3: Three-body recombination event rate K3 (in units 
of cm 6 /s) as a function of the collision energy E (in units 
of mK). The solid lines are the results in Fig. 3 of Ref. |{J. 
The dashed lines are the hyperspherical averages of the R in 
Eq. G3 an d the J = 0, 1, 2 and 3 terms in Eq. lj20l . 



interference between two terms in Eq. In the diago- 
nal term, the last factor peaks at cos 0-3 = \ — \EvjE 
with a width proportional to [(E D /E)(1 + Ed/E)] 1 ' 2 . 
In the limit E Ed, this factor can be approxi- 
mated by a delta function <5(cos 2 a3 — i) multiplied 
by (2tt/V3){E/E d ) 1 / 2 . The resulting expression for R 
agrees with that in Eq. 1(13(1 . The hyperangles a% and a 2 
in Eq. lfT7|) can be expressed as functions of 03 and the 
angle (3% between the vectors ps and p\ — p 2 : 

cos 2 0*1,2 = \ cos 2 a 3 + I sin 2 a 3 

^v 7 ^ cos 0:3 sin 013 cos /33. (19) 



We can get further insights into the factorization ap- 
proximation by considering the angular-momentum de- 
composition of the 3-body recombination rate. The rate 
calculated by Suno et al. in Ref. Q was obtained by 
adding the contributions with total angular momentum 
quantum number J = 0, 1, 2 and 3. These contribu- 
tions are shown in Fig. EI Note that at E = 10 mK, the 
highest energy for which they were calculated, the J = 
and 1 contributions are increasing and they are smaller 
than the J = 2 and 3 terms which are decreasing. The 
angular-momentum decomposition of the rate in Eq. 117(1 
is 
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(20) 



where the Qj(z) are Legendre functions of the second 
kind with branch cuts on the interval —1 < z < +1. 
Their arguments are 

cos 2 a, + \ + E D /E 

Zi = — = . (21) 

V 1 + E D /E cosa^ 

The argument of the Legendre polynomial Pj(z) in 



Eq. is 

COS 2 0-3 — COS 2 CL\ — COS 2 Oi 2 

23,12 = ~ ■ (22) 

2 COS OL\ COS Q-2 

The individual contributions from J = 0, 1, 2, and 3 are 
shown in Fig. [3] In the high energy limit, each individual 
contribution scales like E~ 5 / 2 . Since the total (R) scales 
like £T 2 , the sum over J must provide the additional 
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factor of E 1 / 2 . At E = 10 mK, the J = 0, 1, 2, and 
3 terms are larger than the results of Ref. [t| by factors 
of 711, 12.0, 1.6, and 1.7, respectively. The mismatch 
in the J — 2 and 3 curves in Fig. [3] is small. The large 
mismatch in the J — 1 curves is easy to understand. The 
cross terms in Eq. (1201 are constructive for even J and 
destructive for odd J, and they are much smaller than 
the diagonal terms for J > 3. Thus the J = 1 channel 
is the only one with substantial destructive interference. 
The sum of the diagonal and cross terms is 47% of the 
diagonal terms at 100 mK, 22% at 10 mK, and 4.5% 
at 1 mK. Because of the destructive interference at low 
energies, small corrections to the amplitude can give a 
large correction to the rate. Thus only small corrections 
to the amplitudes for J = 1, 2, and 3 would be required 
for the terms in Eq. (|20|) to match smoothly on to the 
results of Ref. |9|- These corrections presumably come 
from higher orders in the perturbation expansion in the 
exact atom-atom scattering amplitude. 

We now discuss the much larger mismatch between 
the J = curves in Fig. [3] In this channel, there is con- 
structive interference between the diagonal terms and the 
cross terms in the leading order result in Eq. 120fl . The 
corrections from higher orders in the exact atom-atom 
scattering ampitude could bring the curves into agree- 
ment if they interfere destructively with the leading or- 
der term. The destructive interference would have to de- 
crease the J = term at E = 10 mK by a factor of about 
700. We argue that such strong destructive interference is 
plausible by noting that there is strong destructive inter- 
ference in the J = channel in the calculation of Ref. [t| . 
The maximum possible value for K3 at E = can be ob- 
tained from Eq. by setting an = 87.0 A. This value, 
3.7 x 10 -26 cm 6 /s, is indicated by the horizontal arrow in 
Fig. [21 The value of K 3 at E = in Ref. [| is smaller by 
a factor of 52, which indicates strong destructive interfer- 
ence. The J = result in Ref. Q exhibits even stronger 
destructive interference at an energy near 2 mK, where it 
is smaller than at E — by a further factor of 31. Thus 
the suppression from destructive interference at this en- 
ergy is comparable to that required to make the J = 
curves in Fig. [3] match smoothly at 10 mK. A nontriv- 
ial 3-body calculation would be required to verify that 
higher orders in the exact atom-atom scattering ampli- 
tude provide the necessary destructive interference. The 
result would depend not only on the scattering length 
a or ajj, but also on the three-body parameter re* that 
appears in Eqs. @ and |JSJ. 

The production of 4 He dimers, trimers, and tetramers 
has been studied in experiments involving the free expan- 
sion of a jet of cold pressurized 4 He atoms The ex- 
periments are analyzed using coupled rate equations for 
various formation and break-up processes. The rate con- 
stants for processes involving dimers were estimated by 
using interpolations and extrapolations of previous the- 
oretical and experimental results. Our results for dimer 
break-up in Eq. I|10|) and for 3-body recombination in 
Eq. (|15fl can be used to estimate the rates for those pro- 



cesses at large collision energies. The rate constants in 
Refs. [lpj for processes involving trimers were based on 
geometrical estimates only. Our result in Eq. (|12|) for the 
inclusive break-up cross section can be applied to both 
the excited state and the ground state of the trimer. 
Since the binding energy of the ground state trimer is 
about 50 times larger than that of the excited state, the 
condition E ^> Et requires much larger collision energy 
in the case of the ground state trimer. 

Factorization can also be used to calculate the rate for 
3-body recombination of two atoms and a loosely-bound 
dimer into a loosely-bound trimer when the three incom- 
ing particles all have relative momenta much greater than 
h/a. The leading contribution comes from a hard scat- 
tering of the atoms that scatters one of them into a state 
with small momentum relative to the dimer, allowing it 
to subsequently coalesce with the dimer to form a trimer. 
The short-distance factor has a form similar to that in 
Eq. (|13fl . The long-distance factor involves an overlap 
integral of trimer and dimer wavefunctions. The evalu- 
ation of this factor requires a nontrivial 3-body calcula- 
tion, but it is much simpler than the direct calculation of 
the recombination rate by solving the 4-atom Schrodinger 
equation. The factorization approximation could be used 
to calculate the production rate of both the excited state 
and the ground state of the 4 He trimer in the experiment 
of Ref. 

The production of deuterons, 3 He nuclei, and the corre- 
sponding antinuclei have been observed in ultrarelativis- 
tic heavy-ion collisions 0| . These collisions are believed 
to produce a thermalized state with very high energy den- 
sity which, as it expands and cools, makes transitions to 
a quark-gluon plasma and then to a hadron gas. Analogs 
of the factorization formulas in Eqs. i|10|) and (|15|) can be 
applied to the breakup and formation of deuterons and 
antideuterons in the hadron gas phase. 

The strong interactions between atoms with large scat- 
tering lengths implies that high numerical accuracy is re- 
quired to calculate rates for processes involving loosely- 
bound molecules. Using traditional methods, the list of 
processes for which accurate calculations are possible is 
restricted to those involving at most 3 or 4 atoms. By us- 
ing factorization, the list can be expanded to include all 
those for which the long-distance factors involve N < 3 
or 4 atoms. In some cases, such as the break-up cross sec- 
tions in Eq. i|12|) , the long-distance factors can be deter- 
mined analytically and accurate calculations are possible 
even for N > A. 

In the examples of 4 He molecules and the deuteron, 
the large scattering length a arises from a fortuitous fine 
tuning by nature. Another exciting application of fac- 
torization is to alkali atoms near a Feshbach resonance, 
which allows a to be tuned to arbitrarily large values by 
varying a magnetic field. The factorization formulas give 
the leading term in a systematic expansion in powers of 
1 1 a. As the scattering length is tuned to be increasingly 
large, the factorization approximation becomes increas- 
ingly accurate and it applies at increasingly lower colli- 
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